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Abstract 

We present a detailed analysis of the rare exclusive Higgs-boson decays into a sin¬ 
gle vector meson and a photon and investigate the possibility of using these processes 
to probe the light-quark Yukawa couplings. We work with an effective Lagrangian 
with modified Higgs couplings to account for possible new-physics effects in a model- 
independent way. The /i —>■ Hy decay rate is governed by the destructive interference 
of two amplitudes, one of which involves the Higgs coupling to the quark anti-quark 
pair inside the vector meson. We derive this amplitude at next-to-leading order in Og 
using QCD factorization, including the resummation of large logarithmic corrections 
and accounting for the effects of flavor mixing. The high factorization scale g ~ nih 
ensures that our results are rather insensitive to the hadronic parameters character¬ 
izing the light-cone distribution amplitude of the vector meson. The second amplitude 
arises from the loop-induced effective / 177 * and h'yZ* couplings, where the off-shell gauge 
boson converts into the vector meson. We devise a strategy to eliminate theoretical un¬ 
certainties related to this amplitude to almost arbitrary precision. This opens up the 
possibility to probe for 0(1) modifications of the c- and 6 -quark Yukawa couplings and 
0(30) modifications of the s-quark Yukawa coupling in the high-luminosity LHC run. 
In particular, we show that measurements of the ratios Br(/i —?■ T(n5) 7 )/Br(/i —>• 77 ) 
and Br(/i —)• 66 )/Br( 6 , —)• 77 ) can provide complementary information on the real and 
imaginary parts of the 6 -quark Yukawa coupling. More accurate measurements would 
be possible at a future 100 TeV proton-proton collider. 


1 Introduction 


The discovery of the Higgs boson in 2012 [3,3 has established the existence of a new kind 
of elementary particle, which conples to the other particles of the Standard Model (SM) in a 
non-nniversal way. The SM predictions that the Higgs conplings to heavy gange bosons and 
fermions are given by 2m^^/n and ruf/v, where v ~ 246 GeV is the Higgs vacnnm expectation 
valne, have been conhrmed within experimental nncertainties for the W and Z bosons and for 
the third-generation fermions SB- However, no direct measurements of the Higgs couplings 
to the light fermions of the hrst two generations are available at present. It is not difficult to 
come up with models in which these couplings can deviate significantly from those predicted 
in the SM. For example, in B3 it was proposed that the Yukawa couplings may depend in 
a non-trivial way on the Higgs field, and that this might explain the hierarchies seen in the 
spectrum of fermion masses. A more general analysis of different classes of models in which 
the Higgs couplings to fermions can differ significantly from those of the SM was presented 
in 3 . Probing the Higgs couplings to light fermions is thus of paramount importance. This 
includes both flavor-diagonal and flavor-changing interactions. Correlations between the two 
types of couplings, which to some extent are model dependent, have been studied in 3- The 
CMS collaboration has recently reported a slight excess in the search for the flavor-violating 
decay h —)■ 3, which, if interpreted as a signal, corresponds to a branching fraction 

Br(h —)■ = (0.89 %. Not surprisingly, this observation has led to much theoretical 

speculation. 

Our focus in the present work is on the Higgs couplings to light quarks {q ^ t). In a couple 
of beautiful papers, it has recently been proposed that one might get access to these couplings 
by focussing on the rare, exclusive decays h V'y of the Higgs boson [l^, IT], where the 
final state contains a single vector meson V. Such measurements are extremely challenging 
at the LHC, as the corresponding branching fractions are in the range of few times 10“®. 
Nevertheless, observing these processes is not hopeless in view of the fact that in its high- 
luminosity run the LHC will serve as a Higgs factory. With 3 ab“^ of integrated luminosity, 
about 1.7x 10® Higgs bosons per experiment will have been produced [l2| , and even larger Higgs 
samples could be obtained at a future facility such as a 100 TeV proton-proton collider. The 
theoretical description of rare exclusive decays employs the formalism of QCD factorization 


lOMlTj. which was originally developed for the analysis of hard exclusive QCD processes and 


later extended to the more complicated case of exclusive hadronic two-body decays of B 


18L Il9|. In a recent paper, we have systematically developed this approach for the 


mesons 

case of the exclusive decays Z —)■ My, W —)■ My and Z —)■ MW j2^. These processes are less 
susceptible to new physics and can therefore be used to test the QCD factorization approach 
and extract valuable information about the light-cone distribution amplitudes (LCDAs) of 
various mesons M. Observing exclusive radiative decays of heavy electroweak gauge bosons in 
the high-luminosity run of the LHC would provide a proof of principle that this kind of rare- 
decay searches can be performed in a hadron-collider environment. A promising first step in the 
direction of observing the decays Z —)■ J/V' y and Z —)■ Ty, along with the corresponding Higgs 
decays h —)• J/V’y and h Ty, has recently been reported by the ATLAS collaboration |2lj. 

In this work we extend previous studies of exclusive radiative decays of the Higgs boson 
in several important ways. We include QCD radiative corrections and resum large logarithms 
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arising from evolution effects between the Higgs mass scale and a low hadronic scale. We 
comment on the structure of power-suppressed corrections to the factorization formula, study 
the effects of the off-shellness of the photon in the h —)■ 77 * —)• 7 I/ process and include the 
power-suppressed h —?■ 'jZ* —)■ jV contribution. We also take into account the effects of oj—(j) 
mixing. Most importantly, our analysis allows for generic non-standard Higgs couplings to 
SM fermions and gauge bosons, including CP-odd interactions. We devise a strategy which 
allows us to eliminate the theoretical uncertainties related to the dominant h —)■ 77 * —)■ 7 H 
conversion contribution, including possible new-physics effects, to almost arbitrary precision. 
This is a crucial prerequisite for achieving the desired sensitivity to the Yukawa couplings of 
light quarks. Finally, we address all relevant h —)■ H 7 decays to both light and heavy mesons 
in one coherent formalism. For technical details on the QCD factorization approach the reader 
is referred to (ioj ]. 

In the following section we start by dehning effective couplings of the Higgs boson to SM 
fermions and gauge bosons. These include both CP-even and CP-odd couplings. In Section |3] 
we discuss in detail the calculation of the different contributions to the h —)• I /7 decay am¬ 
plitudes in the QCD factorization approach, distinguishing between “direct” contributions 
induced by the Yukawa couplings of the Higgs boson to light quarks (g 7 ^ t) and “indirect 
contributions” resulting from a h —)■ /'yZ* transition followed by the conversion of the 

off-shell boson into a vector meson. We show that uncertainties related to the effective ^77 
coupling strength, as regards to both theoretical uncertainties and possible new-physics con¬ 
tributions, can be eliminated by studying the ratio of the h —)■ I /7 and h —)■ 77 branching 
fractions. In Section H] we present a phenomenological analysis both in the SM and in generic 
new-physics models with modihed Higgs interactions. In particular, we point out that mea¬ 
surements of the ratios of the h T(nS') 7 , h ^ bb and h —?• 77 branching fractions can 
yield highly complementary information on the real and imaginary parts of the bottom-quark 
Yukawa coupling. Section 0 contains a summary of our results and the conclusions. Technical 
details are relegated to hve appendices. 


2 Effective Higgs couplings 

In our analysis we assume SM couplings for all particles other than the Higgs boson. For the 
Higgs couplings to SM quarks and gauge bosons we adopt the effective Lagrangian 
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where sw = sin 6 *vi/ and cw = cos^w are the sine and cosine of the weak mixing angle. We 
use = 0.23126 ± 0.00005 as determined from the neutral-current couplings of the Z boson 


evaluated at /i = rriz | 22 |. = \e^''°‘^Fay is the dual held-strength tensor, and we use a sign 

convention where = 1. Our choice of factoring out a loop factor in the second line is made 
for later convenience. For SM extensions in which the new particles are heavy, the coefficients 
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Figure 1: Two-loo Barr-Zee diagram (left) and effective one-loop contribution (right) 
to the EDM of the electron arising from the CP-odd couplings kf, and k^z in the 
effective Lagrangian ([T]). 


of these higher-dimensional operators are suppressed by two powers of the new-physics scale. 
We emphasize that the above is not a complete list of operators. For instance, we have not 
included higher-dimensional operators of the form and whose coefficients 

are already strongly constrained by data. These operators would enter our analysis only via 
the h —)■ 77 */ 7 Z* one-loop amplitudes, and without loss of generality their effects can be 
absorbed into the coefficients and k^^, k^z- 

Both the CP-even couplings Ki and the CP-odd coefficients ki are real. In the SM kw = 
Kz = Kf = 1, while kf and all the remaining couplings in the second line vanish. Our Kq and 
kq parameters for quarks are related to the corresponding Yukawa couplings by 

-^ = {Kq+ ZKq) — = {Kq+ tKq) — . (2) 


In the last step we have introduced rescaled parameters normalized to the mass of the b quark. 
This will turn out to be a useful definition for the quarks of the first two generations. 

For our analysis we need the quark masses at the scale /i = ruh, where rrih = (125.09 ± 
0.24) GeV is the Higgs-boson mass j^. We define the running quark masses at next-to- 
next-to-leading order (NNLO) in the MS scheme, starting from the low-energy values given 
in j^. This yields rribirnh) = (2.79 ± 0.02) GeV, rricirnh) = (622 ± 12) MeV, ms{mh) = 
(52.8 ± 1.4) MeV, md{mh) = (2.66 ± 0.11) MeV, and mu{mh) = (1.21 ± 0.08) MeV. For the 
top quark we obtain mt{mh) = (166.8 ± 0.7) GeV starting from mtimj) = (163.4 ± 0.7) GeV, 
which we have derived from the present world average obtained in [2^ using the conversion 
tables provided in [^. Present LHG data are largely insensitive to the Yukawa couplings to 


light quark flavors. From a global fit to the measured Higgs rates, the authors of UJ, [26 
have derived the bounds \/|k„P + \ku\‘^ < 1.3 and \/|kd,s,cp + |^d,s,cp < 1-4 at 95% confidence 
level (GL). The corresponding bounds for the original parameters are \/\ku\^ + \ku\‘^ < 3000, 
+ \kd\‘^ < 1500, \/\ks\^ + |ksP < 75 and \/\kc\^ -\- \kc\^ < 6.2. 

Bounds on the GP-violating Higgs couplings to third-generation fermions have been studied 
in jl^. Under the assumption that the Higgs couples to the electron in the standard way 
{K(. = 1, ke = 0), the strongest constraints come from the electric dipole moment (EDM) of 
the electron. They arise from two-loop Barr-Zee diagrams such as the one shown on the left in 
Figure [Hand yield \kt\ < 0.01, |k;,| < 1.9 and \kr\ < 2.4 at 90% GL. Bounds from the neutron 
EDM are approximately one order of magnitude weaker, but they do not reply on assumptions 
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about the Higgs couplings to the electron. The h ^ bb and h —)■ rate measurements at 

the LHC can be used to place upper limits on the combinations which imply 

the stronger bounds (at 95% CL) |Kfe| < 1.44 and |kt-| < 1.24 from CMS j^, and |Kb| < 1.3 
and |kt-| < 1.5 from ATLAS [^. Bounds from the electron EDM can also be derived for the 
local operators multiplied by and R^z in dl])- Evaluating the one-loop contributions to the 
electron EDM shown on the right in Figure [H we obtain in the MS scheme 
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1 - 4s^ ( 3 xz\r].xz\ 


(^7^ T ^e) ■) 


( 3 ) 


where xz = This contribution is logarithmically UV divergent, because the inner 

structure of the effective hyy and h'yZ vertices is not resolved. The term proportional to R^^ He 
agrees with a calculation performed in j^. The subtraction scale jj, should be identified with 
the scale of the new physics, which is responsible for these non-standard interactions. Setting 
/i = Anp = 1 TeV for an estimate, and assuming SM-like Higgs couplings to the electron, 
we obtain from the present experimental bound |de| < 8.7 • 10“^®ecm (at 90% CL) the 
constraint 

\R^^ + 0.09 R^z\ < 0.006 (90% CL). (4) 

Barring a hne tuning of the two contributions, this implies that \R^'^\ < 6 ■ 10“^ and \R'yz\ < 
0.07. If the new-physics scale lies above 1 TeV then these bounds become stronger. With 
Anp = 10 TeV, for example, they improve by a factor of 2. 


3 Radiative hadronic decays of Higgs bosons 


Our focus in this work is on the rare, exclusive radiative decays h ^ V'f, where V denotes a 
vector meson with momentum k and the photon carries momentum q. We will refer to vectors 
orthogonal to the plane spanned by k and q as being transverse. Up to tiny corrections of 
order {rnv/rrihY-, the mass of the vector meson can be set to zero. 

The leading-order Feynman diagrams are shown in Figure |2j In the hrst two graphs, 
the Higgs boson couples to the quark and anti-quark pair from which the meson is formed. 
We refer to this as the “direct” contribution to the decay amplitude 


calculated from first principles using the QCD factorization approach 



It can be 
because the 


energy released to the final-state meson is much larger than the scale of long-distance hadronic 
physics [3,00. At leading power in an expansion in Aqcd/’^/x, the direct contribution can 
be expressed as a convolution of a calculable hard-scattering coefficient with the leading-twist 
LCDA of the vector meson V. The corresponding factorization formula was derived in 0| 
using the formalism of soft-collinear effective theory 0^01 • It was shown that, for a given 
helicity amplitude, the power corrections to the leading term are suppressed by (Aqcd/’^/i)^ 
for light mesons and (rug/m/i)^ for mesons containing heavy quarks of flavor Q. Even for 
the 5-quark these power corrections are negligible. The third diagram in Figure |2] shows a 
different production mechanism, in which the vector meson is produced via the conversion 
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Figure 2: Direct (left and center) and indirect (right) contributions to the h —)• I /7 
decay amplitude. The crossed circle in the third diagram denotes the off-shell h —)■ 77 * 
and h —)■ 7 Z* amplitudes, which in the SM arise first at one-loop order. 


of an off-shell photon or Z boson produced in a h —>■ 77 */ 7 Z* transition Ij. We refer to 
this as the “indirect” contribution. It involves the hadronic matrix element of a local current 
and thus can be expressed in terms of the decay constant /y of the vector meson. The direct 
contribution is sensitive to the Yukawa coupling of the Higgs boson to the quarks which make 
up the vector meson. We shall find that in the SM the direct and indirect contributions to 
the h —)■ 1/7 decay amplitude interfere destructively. They are of similar size for V = T, 
while the direct contributions are smaller than the indirect ones by factors of about 0.06 for 
V = 0.002 for V = cp, and few times 10“^ for V = and u. The sensitivity to the 

Yukawa couplings thus crucially relies on the precision with which the indirect contributions 
can be calculated. We will come back to this point below. 

The most general parametrization of the h —)• Y 7 decay amplitude is 


iA{h —>■ Y 7 ) = 


e/y 


-y 


q-e*yk-e*^ 
k ■ q 




k ■ q 


nV 


( 5 ) 


where both the hnal-state meson and the photon are transversely polarized. From (|5]), the 
decay rate is obtained as 


F(h ^ 1 / 7 ) = 


a/, 


y 


8m/, 


77 + 7 ; 


VI 


( 6 ) 


Here a = 1/137.036 is the hne-structure constant evaluated at = 0 [22j, as appropriate 


for a real photon. We choose to normalize the decay amplitude in ([5]) to the vector-meson 
decay constant /y, which is defined in terms of a matrix element of a local vector current. 
Since we consider neutral, flavor-diagonal mesons, the definition of the decay constants (and 
of other hadronic matrix elements) is complicated by the effects of flavor mixing. In complete 
generality, such a neutral meson V can be regarded as a superposition of flavor states \qq). 
We can thus define flavor-dependent decay constants fy via 


{V {k, e) I |0) = -if^mys*^ ; q = u,d,s,... . 


( 7 ) 


A certain combination of these flavor-specihc decay constants can be measured in the leptonic 
decay V e'^e~. The corresponding decay amplitude involves the matrix element of the 
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(8) 


electromagnetic current Qg We thus define 


Qvfv = Qqfv ) where Qv = ^c]^Qq- 
q q 

Here denote the fiavor coefficients in the naive constituent-quark model, where |p°) = 
[\uu) — \dd)), |a;) = {luu) + \dd)), |0) = |ss) etc. It follows that 

= ^f^ = T.Q',a. = ( 9 ) 

and so on. With these definitions, the electromagnetic decay rate is given by 


r(l/ ^ e+e") = 


4vrg^^/^ 

3mv 


a^{mv) . 


( 10 ) 


In [20|, the so-defined decay constants fv have been extracted using the most recent experi¬ 
mental data. Our motivation for factoring out these decay constants in ([S]) is that they can 
be determined experimentally without adopting any specific model for fiavor mixing. As we 
will discuss in Section |321 the dominant terms in the indirect contributions to the h —?• Hy 
decay amplitudes are proportional to precisely these quantities. 


3.1 Direct contributions to the form factors 


We first consider the direct contribution to the decay amplitude shown in the first two diagrams 
in Figure m in which the Higgs boson couples to the quarks contained inside the vector meson. 
In order to calculate the corresponding contributions to the form factors defined in ([5]), one 
calculates the corresponding partonic amplitudes with on-shell quark and anti-quark states 
and then projects these amplitudes onto the leading-twist LCDA of a transversely polarized 
vector meson V. The relevant light-cone projector reads 38 


My^(/c,X,/i) = + , (11) 

where the dots stand for higher-twist contributions. In analogy with ([7]) and (|H]), we define a 
set of flavor-specific transverse decay constants via 

{V{k,e)\qia>‘-q\a) = (Fe" - Ve*") . (12) 

The quantity fy entering 111 1 1I is then defined by the eombination /^ = (^ Qqfv ) /Qv- 
The transverse decay constants are scale-dependent quantities, since the QCD tensor current 
has a non-zero anomalous dimension. The leading-twist LCDA 0y(a;,p) can be interpreted as 
the amplitude for finding a quark with longitudinal momentum fraction x inside the meson0 
It depends on the choice of the factorization scale /i employed in the factorization formula. 

^Strictly speaking, one should introduce flavor-specific LCDAs and define the product 

fy^{F) F) in terms of a matrix element of a non-local quark current with fiavor q, in analogy with rela¬ 

tion (4) in [20j. Because the LCDAs are normalized to 1, and given the present large uncertainties in the shapes 
of these functions, it is a safe approximation to employ SU(3) symmetry and replace 
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Figure 3: One-loop QCD corrections to the first diagram in Figure |21 Analogous 
corrections exist for the second diagram. 


The two leading-order graphs shown in Figure |2] are supplemented by the diagrams in 
Figure [31 which arise at 0{as)- Including these effects is crucial in order to control the scale 
dependence of the transverse decay constant fy, the Yukawa coupling yq and the LCDA (j)y. 
It will also allow us to resum large logarithms of the form ln(m^//ig)) to all orders in 
perturbation theory. Here /xo ~ 1 GeV is a typical hadronic scale, at which model predictions 
for fy and fy are obtained. We work in dimensional regularization and subtract UV and IR 
divergences in the MS scheme. The product of the bare decay constant times the LCDA of a 
transversely polarized vector meson is related to the product of the corresponding renormalized 
quantities via 

= / dy Zf^{x,y,y) fy{y)(j)y{y,y) , (13) 

Jo 

where at one-loop order 

Z^{x, y, y) = S{x - y) + V^{x, y) + 0{al) , (14) 

with Cp = 4/3. The relevant one-loop Brodsky-Lepage kernel reads BQ 




^ , 0{.x-y) 

a^(l - y) -+ 2/(1 - x) 


y-x 


x-y 


(15) 


For the decays h —?• Ry, which are mediated by (pseudo-)scalar currents, an overall UV 
divergence remains, which is cancelled by the counterterm for the Yukawa coupling, derived 
from 2/g,bare = y^Zy{fi) yq{y) with 

Zy{p) = I + 0(al). (16) 

When dealing with pseudo-scalar currents we employ the’t Hooft-Veltman (HV) scheme [doj ]. 
in which 75 = anti-commutes with the four matrices y^ with /i G {0,1, 2, 3}, while it 

commutes with the remaining (d—4) Dirac matrices y(|(. While this definition is mathematically 
consistent, it violates the Ward identities of chiral gaug e theories by finite terms, which must 
be restored order by order in perturbation theory This is accomplished by performing 

the finite renormalization P = Z^yPhy of the pseudo-scalar current P = q'y^q, where 

= 1 - 2 Cf Yd + 0(ai). (17) 

TT 
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By evaluating the relevant Feynman graphs in Figures [2] and [31 we hnd that the direct 
contributions to the form factors in the amplitude decomposition (|5]) are given by 


-indirect - I^V QvFv , -^ 24 irect “ QvFy , 


V 


where we have defined 


- _ _J_ ^ fv^ 

KV — 2_^ l^qQq 


1 £± ’ 




KV — 
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— V ^ Q 


fv 


(18) 


(19) 


The reduced form factors Fy are given by 
mbifJ-) /y(/i) 


Fy = 


fv 




'0 


x(l-x') I ^ ^ ^ 


( 20 ) 


with ^ 

h{x, mh, /u) = 2 In [(r(l — x)] ^In + In^ x + ln^(l — x) — 3 . (21) 

This result agrees with a previous calculation performed in apart from a typoH We 
have expressed the Yukawa coupling in terms of the running 6-quark mass using the second 
relation in (|2]). We focus primarily on the cases V = J/ip and T{nS), where to an excellent 
approximation the vector meson contains a single quark flavor q, and hence kj/^ ~ Kc and 
kT(nS) ~ i^b, and similarly for the CP-odd parameters ky. For the light mesons V = p^, u and 
(j), on the other hand, flavor-mixing effects can be important. This concerns, in particular, the 
possibility of a small admixture of an |ss) flavor component in the wave functions of and 
u, which can be important due to the smallness of the Yukawa couplings to the up and down 
quarks in the SM. Since the meson is a member of an isospin triplet, its flavor mixing into 
|ss) can only be caused by electromagnetic interactions or isospin-violating effects in QCD. 
Both types of effects are estimated to be very small, and hence we expect that \fpf/f^\ 1- 
To good approximation we can thus use the naive relation 


KpO ^ 


2 ^-^ bid 


SM 

-)■ 


6.1 ■ 10 "^ 


( 22 ) 


The situation is different for the case of the uj meson, whose mixing into an |ss) flavor state 
can be non-negligible. In Appendix [^ we derive explicit expressions for the parameters k^^ 
and k^ in a simple flavor-mixing scheme for the u—(j) system. Assuming that |ks| ^ \ku^d\ like 
in the SM, and working in the SU(3) limit and to hrst order in the small mixing angle 0^^, 
we obtain 


Ku 


K<p K, Ks 


2ku-kdFV2kse^^{ml) ^ (-0.08 + 26.8 0 ,^)-10-3, 


1 + 


x /2 


SM 


^ 0.019 + 0.013 0,,^. 


(23) 


^These authors use the pole mass instead of the running quark mass in the prefactor, which adds 
—31n(/x^/m^) — 4 to the kernel/i(a:, m/j,/r). In eq. (130) of one finds instead—3ln[/Li^/(—m^)] — 4. We are 
grateful to the authors for confirming this mistake. 















In the SM the contributions from the up and down quarks almost precisely cancel in and 
hence the contribution induced by cn —0 mixing is likely to be the dominant one. Existing 
estimates for the mixing angle 9^^ derived from mass-independent analyses include 9^^ 0.05 


42| and 9^^ ^ 0.06 [43|, [44 


On the other hand, in a more recent mass-dependent analysis 
the values ~ 0.008 and 9i^(f,{m^) ~ 0.081 were obtained [i^. We conclude from this 

discussion that ~ Rg to good approximation, while a more accurate description of flavor¬ 
mixing effects would be required before the quantity R^j can be interpreted reliably in terms 
of quark Yukawa couplings. 

In the factorization formula fl20ll all non-perturbative hadronic physics is contained in the 
decay constants and the LCDA. The quantity multiplying the LCDA under the integral is 
the hard-scattering coefficient, which can be calculated in perturbation theory. It depends 
on the momentum distribution of the quark inside the hadronic bound state. QCD-based 
model calculations of LCDAs are typically performed at a low hadronic scale /xo ~ 1 GeV. If 
such a low value is chosen for the factorization scale /r in fl2n|l . the hard-scattering coefficient 
contains large logarithms of the form (og ln(m^//iQ))” with \n(m1 /~ 9.7, which should be 
resummed to all orders in perturbation theory. To perform this resummation, it is convenient 


to use the expansion of the LCDA in the basis of Gegenbauer polynomials, which reads 13|, [17 


4>y{x, fi) = 6x(l — x) 




n=l 


(24) 


We can then rewrite fl20|) as a sum over Gegenbauer moments, finding 


Fv = 


Qmbifx) /y (/r) 
V fv 


1 - 




71 




hi'mi,) , 


(25) 


with 


A'(mft) = ti) 


(26) 


n=0 


The factor in brackets in fl25|l precisely compensates the scale dependence of the product 
Tjihi u) /y ()w), while the quantity is formally scale invariant. Using a technique developed in 
(20[, we obtain for the hard-scattering coefficients in moment space the closed expression 


Cn{mh,n) = 1 -h 


Att 


-4(i7„+i-l) ( 


ml 


In ^ — fvr I AHf 
\ M 


n+1 


3 -|- AiTT 


+ 0{al), (27) 


where Hn+i = Xlfcii i ^^e harmonic numbers. As a consequence of the symmetry of the 
hard-scattering coefficient under the exchange x -f-)- (1 — x), the sum in fl26ll runs over even 

a^ln-^) can now be resummed 
readily by choosing the factorization scale of order /i ~ mh and evolving the scale-dependent 
quantities /y(/i) and up to that scale. The solution of the corresponding 

renormalization-group (RG) equations at NLO of RG-improved perturbation theory is dis¬ 
cussed in Appendix [HI At leading order in QGD the Gegenbauer moments a)(^(/i) in fl2Tj) are 
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Meson V 

fv [MeV] 

fi,{2GeV)/fv 

ay^(/io) 

Qv 

Vy 


216.3 ± 1.3 

0.72 ± 0.04 

0.14 ±0.06 

1 

F2 

^w) 

u 

194.2 ±2.1 

0.71 ±0.05 

0.15 ±0.07 

1 

3V2 

_ 

3\/2 

0 

223.0 ± 1.4 

0.76 ± 0.04 

0.14 ±0.07 

1 

3 

1 <52 

4 3 


Table 1: Hadronic input parameters for light vector mesons. The decay constants fv 
are extracted from data on the electromagnetic decay widths V —)■ jl^, while the 
ratios fy/fv are derived from a compilation of theoretical predictions. The values of 


49 . The last 


the Gegenbauer moments at the scale fio = 1 GeV are taken from [ 
two columns show the effective charges Qy and vy dehned in ([ 8 ]) and below fl32ll . 


renormalized multiplicatively, such that 


13, 30 


^ , with = 8Cf (i^n+i - 1) • (28) 

V Ois (/Tq) / 

Here /iq denotes a low reference scale, while /x = 0{mh) is the hard factorization scale, to which 
the LGDAs are evolved. The NLO corrections to these relations are discussed in Appendix iBl 
They have a negligible impact on our numerical results. Ah of the anomalous dimensions are 
strictly positive (for n 7 ^ 0), which implies that a)(^(/i) —?• 0 in the formal limit /x —?• cxo. In this 
limit the leading-twist LGDAs approach the asymptotic form 6 a:(l — x). Similarly, it follows 
from relation flB. 2 D in Appendix [B] that the transverse decay constants of vector mesons vanish 
in the asymptotic limit, i.e. /y (^— )■ 0 for /x —)■ 00 . 

It has been emphasized in [20(| that RG evolution effects render our predictions rather 
insensitive to the precise values of the Gegenbauer moments. From fl26|) . we obtain 

Re/y(m/i) = 1.01 -|- 1.13ay^(mft) -|- 1 . 21 ap^(m/i) -|- 1.29a^-^(mft) -|- 1.35a^-^(m/i) -|- ... 

(29) 

^ 1.01 + 0.51a^^(/io) + 0.36a]f^(/io) + 0.29a^^(/io) + 0.24a^^(/io) + ... . 


While ah Gegenbauer moments have 0{1) coefficients at the high-energy scale, the coefficients 
of the higher moments are strongly reduced when one expresses the answer in terms of moments 
normalized at the low scale /xq = 1 GeV. 

In order to obtain numerical predictions for the reduced form factors we need as hadronic 
input parameters the decay constants fy and fy and the Gegenbauer moments of the 
various vector mesons. As mentioned earlier, the decay constants fy can be extracted from 
experimental data, and up-to-date values have been derived in . The ratios fy / fy needed 
in fl25l) must be obtained using some non-perturbative approach, such as lattice QGD, light- 
con e Q GD sum rules or the non-relativistic effective theory NRQGD for heavy quarkonia 
46l. l47l | , which provides a systematic expansion of hadronic matrix elements in powers of the 
small velocity v ~ as{mQv) of the heavy quark in the quarkonium rest frame. Details of 
such determinations are reviewed in Appendix O In Table [1] we compile the relevant input 
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Meson V 

fv [MeV] 

/j^(2GeV)//,,- 

CTy(/io) 

Qv 

Vy 

J/^ 

403.3 ±5.1 

0.91 ±0.14 

0.228 ±0.005 ±0.057 

2 

3 

1 

1 

T(l^) 

684.4 ± 4.6 

1.09 ±0.04 

0.112 ±0.004 ±0.028 

1 

3 

1 1 

4 3 

T(2S) 

475.8 ±4.3 

1.08 ±0.05 

0.144 ±0.007 ±0.036 

1 

3 

1 1 

4 3 

T(3S) 

411.3 ±3.7 

1.07 ±0.05 

0.162 ±0.010 ±0.041 

1 

3 

-1 o2 

1 _|_ 

4 3 


Table 2: Hadronic input parameters for heavy quarkonium states. The decay constants 
fv are extracted from data on the electromagnetic decay widths V —?• l~^l~ j^, while 
the ratios fy/fv are derived from NRQCD scaling relations. The width parameters 
ay are obtained from relation fl3Up . where the first error is of parametric origin and 
the second one parameterizes the uncertainty due to higher-order effects. The last two 
columns show the effective charges Qy and vy. 


parameters for light vector mesons. Because of the lack of information about higher Gegen- 
bauer moments we can only keep few terms in the infinite sum ([23]). The systematics of the 
Gegenbauer expansion has been discussed in [^, where it was pointed out that the higher 
moments with n S> 1 fall off faster than 1/n. Indeed, high-rank Gegenbauer polynomials 
Cn^'^\‘2x — 1) with 77, ^ 1 would resolve structures on scales Ax ~ 1/n. For a light vec¬ 
tor meson V, it is reasonable to assume that the LGDA 4>y{x) does not exhibit pronounced 
structures at scales much smaller than 0{1). To estimate the impact of higher moments we 
use ap^(po) = ±0.15 for our error estimates. Relation fl2^ suggests that the effect of yet 
higher-order terms is small. 

The LGDAs of heavy mesons exhibit a different behavior, since the presence of the heavy- 
quark mass introduces a new scale. For a quarkonium state V ~ (QQ) composed of two 
identical heavy quarks, the LGDA peaks at x = 1/2 and has a width that tends to zero in 
the limit of inhnite heavy-quark mass. The second moment of the LGDA around x = 1/2 can 
be related to a local NRQGD matrix element called {v‘^)v Including the one-loop QGD 
corrections calculated in [3^, we obtain 



if 


{v‘^)v Cpasifi) I'28 2 

3 dvr V9 3 J 


(30) 


A critical discussion of the extraction of the parameters {vfv for different quarkonium states is 
presented in Appendix O Using the values compiled there, but with increased error estimates, 
we obtain the ratios of decay constants and the width parameters cry(juo) at the low scale 


juo = 1 GeV shown in Table [2j As a reasonable model at the scale /uq we adopt the form 20 


0y(x,/7o) = N, 


4x(l — x) 

\p2^Gy 


exp 


1\2 


(^- 2 ) 

2a^ 


(31) 


where the polynomial in front of the Gaussian factor ensures that the LGDA vanishes at the 
endpoints x = 0,1. In order to estimate the uncertainties related to the functional form and to 
capture the effects of unknown higher-order corrections to relation (I5U]1 . we include a second 
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Meson 

Form factor with errors [%] 

Gombined value [%] 

Fpo 

4.30 ; p ± 0.03^, ± 0.24/ ± 0.12,, ± 0.22,, 
±*(0.67l°;J^^ ± 0.00^, ± 0.04/ ± 0.03,, ± 0.06„) 

(4.30 ±0.35) ±*(0.67 ±0.14) 

iG 

4.26 M ± 0.03^, ± 0.30/ ± 0.14„ ± 0.21„ 

±*(0.66± 0.00^, ± 0.05/ ± 0.03,, ± 0.06„) 

(4.26 ±0.40) ±*(0.66 ±0.14) 

F<p 

4.53 M ± 0.03^, ± 0.24/ ± 0.15„ ± 0.23„ 

±*(0.70 1 p ± 0.01^, ± 0.04/ ± 0.04,, ± 0.06„) 

(4.53 ±0.37) ±*(0.70 ±0.15) 

Fjjp, 

4.54 ^± 0.03^, ± 0.70/ 

±*(0.63 p ± 0 .00^, ± 0 .10/1 

(4.54 ± 0.72) ±i(0.63± 0.14) 

Ft{\s) 

5.o4;“”„±o.o4„.±o.i8,i!>:»?„ 

+i( 0.66 1 „ ± 0 , 00 „. ± 0 . 02 ; ,v) 

(5.04 ± 0.21) ±i(0.66± 0.10) 

Fr{2S) 

5.09 „ ± 0.04„. ± 0.24, 1 „ 

+i(0.68 „ ± 0.00„. ± 0.03, 

(5.09 ±0.27) ±*(0.68 ±0.11) 

Ft{3S) 

5.ii1S:SF±')-')4„.±o.24,1S:15„ 

+i(0.69 „ ± 0.00„. ± 0.03, 

(5.11 ±0.29) ±*(0.69 ±0.12) 


Table 3: Theory predictions for the reduced form factors Fv including error estimates. 


error of ±25% on the ay parameters. Given this form, we compute the hrst 20 Gegenbauer 
moments at the low scale evolve them up to the factorization scale /i ~ mh using fl28|) . 
and use these results in evaluating relation fl 2 B]) . 

We are now ready to present our numerical results for the direct contributions to the 
reduced form factors Fy in fl25|l for various vector mesons, including detailed error estimates. 
They are collected in Table |3l The different sources of theoretical errors contain a perturbative 
uncertainty (subscript “/i”), which we determine by varying the factorization scale /i between 
m/j /2 and 2 m/i. Once the NLO corrections are included our results are very stable under 
scale variations. The scale uncertainties are larger for the imaginary parts than for the real 
parts of the form factors, since these start at 0{as) and there is thus no compensation of 
the scale dependence. We emphasize, however, that the imaginary parts only have a small 
impact on our numerical predictions for the decay rates. We also include the uncertainty in 
the value of the 6 -quark mass, which has a very small impact. The uncertainties related to 
hadronic parameters include the ratio fy / fy (subscript “/”) and uncertainties in the shapes 
of the LGDAs, as modeled by the values of the Gegenbauer moments and ap- (light 
mesons) and the width parameter ay (heavy mesons). These hadronic uncertainties are the 
dominant sources of errors. The last column in the table shows the results obtained when all 
errors are added in quadrature. These numbers will be used for our phenomenological analysis 
in Section 01 We observe that the spread of the results for the form factors Fy for different 
vector mesons is rather small. The theoretical uncertainties on the real part of Fy are typically 
between 4% and 9%. The only exception is Fj/^p, for which the uncertainty in the ratio of 
decay constants is about 16%. It would probably be possible to reduce this uncertainty by 
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Figure 4: Examples of electroweak radiative corrections (top row) and higher-order 
QCD radiative corrections (bottom row) to the h —)■ E 7 decay amplitudes. 


performing a more detailed NRQCD analysis. 

Power-suppressed corrections to our results (120]) and fl2S]) can be org anized in an expansion 
in (Aqcd/''^/i)^ for light mesons and {mv/rnhY for heavy mesons [20|. They are at most of 
order 10“^ for V = T(nS), 10“^ for V = J/V’, and 10“^^ for all light mesons. It is thus safe 
to work with the leading-order terms. In our analysis we neglect two-loop QCD corrections, 
whose effects should be covered by the error we estimate from scale variations, and one-loop 
QED or electroweak radiative corrections, a few examples of which are shown in the top 
row in Figure 01 For flavor-diagonal final-state mesons, the first diagram involves a factor 
^ ~ 0.7 ■ 10“^ instead of the Yukawa coupling yg in the diagrams in Figure [2], while the 
second diagram involves a factor ^ ~ 2 ■ 10“^. This is smaller (by roughly a factor 10) 
than the charm-quark Yukawa coupling and of the same order as the strange-quark Yukawa 
coupling. If the goal is to reach sensitivity to the strange-quark Yukawa coupling in the SM, 
then these electroweak corrections should be calculated. However, such a level of sensitivity 
will be out of reach at the LHC. In the bottom row in Figure 0] we show other examples of 
neglected diagrams, which involve ah—?- gg{'y) transition followed by the conversion of the 
two gluons into the final-state meson. The diagram on the left corresponds to a two-gluon 
LCDA of the vector meson, which does not exist at leading twist due to the Landau-Yang 
theorem. The graph on the right is analogous to the second diagram shown in the first row, 
but with the internal W bosons replaced by gluons. A naive estimate indicates that the two 
types of effects should be of similar magnitude. 

3.2 Indirect contributions to the form factors 

We now proceed to study the photon- and Z-pole contributions to the decay amplitude shown 
in the third diagram in Figure [2j The crossed circle in this diagram represents the off-shell 
h —?• 77 * and h —?• 'yZ* amplitudes, which in the SM are induced by loop graphs involving 
a virtual charged fermion or a VF boson (in unitary gauge), as shown in Figured Since the 
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Cx] 


Figure 5: One-loop SM contributions to the effective h 77 and h'^Z vertices. 


indirect contributions to the decay amplitudes are numerically dominant over the direct ones 
and our goal is to compute them with the highest possible accuracy, we include the effect that 
the intermediate gauge boson is slightly off shell (/c^ = ruy), and we keep the full dependence 
on the meson mass even though this is a very small effect. 

The exact one-loop expressions for the off-shell h —)■ 77 * and h —)■ 7 Z* amplitudes have 
been derived in jlH. Using these results and extending them to the case of the pseudo-scalar 
Higgs couplings in ([T]), we find 


1,indirect 

2,indirect 


aijnv) — my 

71 my V 

. a{mv) m\ — my 

% 

TT my V 


Qv C^^iyXy) 
Qy C^^yy) 


Vy my 

[Sy/Cw)^ ~ 

Vy my 

{syrCw)^ ~ 


C^z{xy) , 
C^zixy) , 


(32) 


where xy = my/m\ accounts for the effects of the off-shell boson, and vy = is defined 

in analogy with Qy in ([ 8 ]), where Vf = — SyrQf are the vector couplings of the Z boson 

to fermions. It is a safe approximation to neglect flavor-mixing effects for the subleading 
contribution from h —)■ 7 Z* —)■ 7 U. At one-loop order, the loop functions are given by 


Cj^{xy) 

C^z{xy) 

and 


X] S -^ Xy) + ^ AfQi, Xy) - ^ Xy) + , 


E 


Ka 


2NcQgVq 


Af{Tq, Xy) + ^ Af{Ti,Xy) - ^ {tw, Xy) + K^z , 


KW 

2 


C'y'yi^Xy ) ^ ^ ITiqN qQ q Xy ) -|- ^ ^ Kl Qi BjQi, Xy ) -|- Kjy , 

q I 

C-,z(xv) = ^ ii,NcQ,v, Bi(t„ Xv) + Si Qm Bi(ti,Xv) + . 


(33) 


(34) 


The first two terms in each coefficient are the contributions from the quarks and leptons, 
the third term in and C-^z arises from gauge-boson loops, and the last term accounts for 
possible new-physics contributions parameterized by the operators shown in the second line 
of ([I]). We have introduced the dimensionless variables Tf = Am'j/m'l (for f = q,l) and tw = 
4myx/m\. We use the running quark masses mq{mh) when evaluating the variables Tq, which 
is appropriate in view of the large momentum transfer in the loop. Explicit expressions for the 
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loop functions Af, A'^ and Bf are given in Appendix [Dl In the SM we have Kg = ki = kw = 1 
and = K^z = 0. The effective Higgs couplings ki entering in fl3T|l all vanish in the SM. 

Since for small values of r/ the fermion loop functions Af{Tf,xv) and Bf{Tf,xv) are pro¬ 
portional to Tf, it suffices for all practical purposes to keep the contributions from the third- 
generation fermions. The effects of the off-shellness of the photon that converts into the final- 
state vector meson gives rise to very small corrections. The relevant variable xy = my/ml 
varies between 3.8 x 10“® for V = p and 6.9 x 10“^ for T(3S'). Note also that the contribution 
of the h —)■ ■yZ* amplitude in fl3^ is by itself strongly power suppressed. Numerically, we 
obtain 

^ 77 ( 0 ) = ^77 “ 4:A64kw + 0.920^4 — (0.012 — O.Ollf)^,- — (0.007 — O.OOSi)^;, 

- (0.015 - 0.010i)kc - 0.001k, -f ... , (35) 

C'^z(O) = K^z - 2A73 kw + 0.132Ki - (0.004 - O.OOdi)^;, - (0.002 - 0.00H)kc + ... , 


and similar expressions hold for the CP-odd coefficients C-yy and Cyz- Notice that the con¬ 
tributions from second-generation fermions are very small even if one assumes that their 
Yukawa couplings are as large as those of the Yquark (i.e., for kc,, = 1). In the SM, we have 
^ 77 ( 0 ) = —3.266 -|- 0.021i and CyziS^) = —2.046 -|- 0.005f, while the CP-odd coefficients C. 


77 

and Cyz vanish. In the expressions for the form factors in 0321) these coefficients are weighted 
by different factors. Note also that to good approximation the indirect contributions to the 
form factors are proportional to 1 /my and hence they are larger for lighter mesons. 

Let us now discuss the impact of QCD and electroweak radiative corrections on the above 
results. Gluon exchanges between the two quark lines that make up the final-state meson in 
the third diagram in Figure [2] are of non-perturbative nature and are accounted for by the 
meson decay constant /y. Radiative corrections connecting the effective hyy and hyZ vertices 
in the last diagram in Figure [2] to the final-state quarks give rise to graphs which no longer 
receive the enhancement proportional to m\/my from the photon propagator. Also, at least 
two gluons would need to be exchanged by color conservation. The effect of such diagrams is 
negligible. We thus only need to worry about QCD corrections to the h —>■ 77 * and h —)■ yZ* 
amplitudes, which arise from two-loop graphs in which a gluon is exchanged inside the quark 
loop. These corrections have been calculated in numerical form in 5^ and analytically in 
In practice the corrections are only relevant for the top-quark contribution. Their effect 

Two-loop electroweak corrections to 
(see also |^). One finds that 


loop. 

13 . 


is to enhance the decay amplitude by a few percent, 
the k —)■ 77 amplitude in the SM were calculated in 5^, 1551 
they are small and negative. For rrih = 125.09 GeV, the effects of QGD and electroweak 
corrections nearly cancel each other, leaving a total correction of about —0.2% In our 
phenomenological analysis we will ignore radiative corrections on the indirect contributions 
computed here, as their combined effect is well below the 1% level. We will device a strategy 
where the dominant contribution proportional to C'.y.^(O), including radiative corrections, drops 
out. Radiative corrections then only have a tiny impact of the coefficients of the subleading 
terms. 
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3.3 Reduction of theoretical uncertainties 


If we assume that the Higgs couplings to the electron are close to those predicted by the SM, 
the CP-odd form factors are more than two orders of magnitude smaller than . For 
a hrst discussion we can thus focus on the form factors F^■ Keeping only the numerically 
signihcant terms, we hnd 

= 0 . 022 kw - O.OOSKi - O.OOSk^^ - (0.017 ± 0.001)^6 + ... , 

= -O.mnw + 0.030^4 + 0.033^.^.^ + (0.030 ± 0.005)hc + ... , (36) 

Ff = 0.206kvf - 0M5Kt - 0.049^.^.^ - (0.015 ± 0.001)k^ + ... , 


where the last term in each line represents the direct contribution. We have dropped the 
small imaginary parts of the latter, whose impact is tiny. Replacing ^ Kg ^ O.OIOks and 
Kc ~ 0.223fi:c one obtains equivalent expressions in which the modihed Higgs couplings are 
expressed as corrections to the SM Yukawa couplings. The challenge is to detect the small 
impact of the direct contributions in the last two cases. 

To this end, it is essential to have absolute conhdence in the precision with which the indi¬ 
rect contributions can be calculated in the SM, and to be able to subtract these contributions 
in a reliable way without assuming that the SM is correct. The latter task can be accomplished 
because the off-shellness of the photon in the h —?• 77 * contribution as well as the h —?• 'yZ* 
contribution in the third graph in Figure [2] are both very small effects. It is therefore possible 
to eliminate the main dependence of the indirect contributions on the new-physics parameters 
by considering the following ratio of decay rates: 

Br(h -)■ V'y) _ T{h -)■ I/ 7 ) _ 87ra‘^{mv) Qvfv A _ V + I’^cp - ^v\‘^ 

Br(h ^ 77 ) F(h—>- 77 ) a rriy \ ml J 1-|-|rcpp 

(37) 

Taking such a ratio has the additional advantage that one becomes insensitive to the unknown 
total width of the Higgs boson, and hence one obtains directly the ratio of branching fractions. 
One can even go one step further and eliminate the sensitivity to the decay constant fv by 
using dTnjl and considering the ratio 

mv Br(h ^ 1 / 7 ) _ 6 / _ V |l - Ay|^-F |rcp - 
F(l/^ e+e“) Br(h —)■ 77 ) a\ ml) l-|-|rcpP 

The only remaining hadronic uncertainties are now contained in the calculation of the reduced 
form factors Fv, which we have collected in Table |3l 

The explicit expressions for the various quantities entering the right-hand side of 0371) are 
rcp = C^^{0)/C^^{0) and 


Ay — —Kv 


Fv nmyv 
a{mv)rn\ 


Ay — —Ky 


Fy TTmyV 
C'.^.y(O) a{my)m\ 


C^^{xv) - 077 ( 0 ) ^ vy C^z{^) 

<^ 77 ( 0 ) <^ 77 ( 0 ) ’ 

C^^{xv) - (F^^(O) rr^ vy C^z{^) 

C'.y.^(O) m? 2 , Qvi^w^w)'^ ^77(0) 


(39) 
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where we work to leading order in the small ratios and xy = m\r/m\. Since the 

individual terms in these expressions are all normalized to C'^^(O), it is convenient to dehne an 
effective parameter k®® by normalizing this coefficient to its SM value. Specifically, we write 




1.275KH' - 0.282k, + (0.004 - 0.003!)k^ + (0.002 - 0.002i)Kt 


+ (0.004 — 0.003i)fi:c 


0.306k. 


77 


/(I 


0.0060 • 


(40) 


Here and below the omitted contributions have coefficients equal to zero up to the indicated 
number of digits. From the empirical fact that the Higgs couplings to W bosons and to the 
third-generation fermions agree with their SM values within errors, it follows that k®® cannot 
be too different from its SM value, except perhaps for a possible new-physics contribution 
parameterized by Note that a tiny imaginary part of k®®, which can be present for 

non-standard values of the k* parameters, would have no noticeable impact on our analysis. 
The ratio rcp entering in flTTIl vanishes in the SM and is entirely due to the various CP-odd 
new-physics parameters kj. We find 


rev = 


-(0.429 + 0.003i)kt + (0.004 - 0.003i)k^ + (0.002 - 0.002i)kb 


k: 


eff 

77 


(0.005 - 0.003i)kc - (0.306 + 0.002i)k, 


(41) 


77 


k: 


eff 

77 


Under the assumption that the Higgs couplings to the electron are approximately SM like, 
the upper bounds |kf| < 0.01 and < 0.006 mentioned above relation ([3]) imply that the 
first and last term in the numerator, which have the largest coefficients by far, are at most 
0.004 and 0.002 in magnitude, respectively. It follows that the magnitude of rcp can at most 
be of order 1%, and hence the impact of this parameter in fl3^ is likely to be negligible. 
We emphasize that the parameter rcp can in principle be probed experimentally by studying 
h —)■ 77 decays, in which both photons undergo nuclear conversion to electron-positron pairs 
3 . In practice, such a measurement appears to be very challenging. 

We now present our numerical results for the CP-even coefficients Ay for the various 
mesons. The complete expressions are collected in Appendix [El They contain direct contri¬ 
butions proportional to the relevant kq parameters and indirect contributions, which are due 
to the power-suppressed h —)■ 7 Z* —)■ 7 U contribution and the effect of the off-shellness of 
the photon in the h —)■ 77 * —)■ 7 U contribution. These latter terms are significantly smaller 
than the theoretical uncertainties in the direct terms. In our phenomenological analysis we 
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will keep these small effects but evaluate them in the SM. This gives rise to the expressions 


Apo = 


A. = 
A^ = 
Aj/^ = 


(0.068 ± 0.006) + i( 0.011 ± 0 . 002 ) 
(0.068 ± 0.006) + i( 0.011 ± 0 . 002 ) 
(0.093 ± 0.008) + i(0.015 ± 0.003) 
(0.281 ± 0.045) + i(0.040 ± 0.009) 


KpO 


K 


eff 


77 


k: 


eff 


77 


j k: 


eff 

77 


+ 0 . 00002 , 


0 . 00011 , 


+ 0.00014, 


(42) 


Kr 


K] 


eff 

77 


+ 0.00005, 


and 


A 


T{1S) 


(0.948 ± 0.040) + i(0.130 ± 0.019) 


-^ +0.0184- 0.0015i, 


K] 


eff 

77 


Ax(25) = [(1.014 ± 0.054) + i(0.141 ± 0 . 022 )] + 0.0207 - 0.0015i, 


(43) 


A 


T(3S) 


(1.052 ± 0.060) + i(0.148 ± 0.025) 


77 


eff 


+ 0.0221 - 0.0015+ 


K. 


77 


Approximate expressions for RpO, and R^ have been given in fl22D and fl2^ . The constant 
terms in the above results show the tiny power-suppressed corrections. Only for the T(nS') 
states they reach the level of percent. Our complete expressions for the CP-odd coefficients Ay 
are also given in Appendix [El It is a good approximation to only keep the direct contributions 
in these terms, which are likely to give rise to the dominant effects. Their coefficients are the 
same as in the expressions above, but with Rg replaced by Rg and Kb replaced by kb¬ 
it is interesting to compare our result for the quantities Ay with corresponding expressions 


obtained by other authors. From [ll[ one can extract ApO = (0.095 ± 0.020) {2ku + Kd)/S, 
A^ = (0.092 ± 0.021) {2Ru + Rd) and A^ = (0.130 ± 0.027)Rs, while from one can obtain 
Aj/-fp = (0.392i0.053)Kc; Ax(i 5 ) = (1.048 + 0.046)/i)6, Ax( 25 ) — (1.138 + 0.053)Kf, and Ax( 35 ) = 
(1.175 + 0.056)Kf,. These values are systematically higher than ours due to the fact that these 
authors have not (or not fully) included QCD radiative corrections and RG evolution effects 
in the direct contributions. For the T(nS') states it is important to keep the small imaginary 
parts of the direct contributions, since in the SM the real parts almost perfectly cancel in the 


combinations 1 —Ay in fl37jl . The result for A^^ obtained in Uj misses the contribution from 


UJ 


(f) mixing and contains a sign mistake in front of Kd- Note also that our predictions for the 
Ay parameters of light mesons are significantly more accurate than those obtained in (llj |. 


4 Phenomenological results 

We begin by quoting our benchmark results for the h —)■ Ry branching fractions in the SM. 
For a Higgs mass of ruh = (125.09 + 0.024) GeV, the SM value of the h —)• 77 branching ratio 
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is (2.28 ± 0.11) ■ 10 ^ Using this result, we obtain for the decays into light vector mesons 


Br{h p^-f) = ( 1.68 ± 0 . 02 /^ ± 0 . 08 /^^.^.^) ■ 10 "® , 

Br(h ^ u-f) = (1.48 ± 0.03/^ ± 0.07,,^.^.^) ■ 10"® 

Br(h ^ 07) = (2.31 ± 0.03/^ ± 0.11/,^.^.^) ■ 10"® , 


(44) 


where we quote separately the uncertainties due to the vector-meson decay constant /y and the 
h —)■ 77 branching ratio, the latter being the dominant source of uncertainty. Our predictions 


are systematically lower and more accurate than those obtained in [llj], where the values 
Br(h —)■ p° 7 ) = ( 1.9 ± 0 . 15 ) • 10 “®, Br(h —)■ u'y) = ( 1.6 ± 0 . 17 ) • 10 “® and Br(h —)■ 07) = 
( 3.0 ± 0 . 13 ) ■ 10 “® are quoted. While the first two results are compatible with ours within 
errors, there is a significant difference for the important mode h —)■ 07. For decays into heavy 
vector mesons, we find 


Br(h ^ J/07) = (2.95 ± 0.07^,/^ ± 0.06dyect ± 0.14;,^.,.,) ■ 10"® , 

Br(h ^ T(l^) 7 ) = (4.61 ± 0.06/^(,,j 1 li? direct ± 0.22;,^^^) ■ lO"" , 

Br(h ^ T(2^) 7 ) = (2.34 ± 0.04;^^^^^ direct ± 0.11,,^.,.,) ■ 10"^ , 

Br(h ^ T(35) 7) = (2.13 ± 0.04y^(3^^ direct ± 0.10;,^.,.,) ■ lO’^ . 


In these cases there is an extra source of theoretical uncertainty related to the calculation of the 
direct contribution to the decay amplitude. Note that there is an almost perfect cancellation 
between the direct and indirect contributions to the h —)■ T {nS) 7 decay amplitudes, and as 
a consequence the resulting branching ratios are roughly three orders of magnitude smaller 
than the h —)■ J/07 branching fraction. For comparison, we note that the branching ratios 
found in [s^ read (2.79 ■ 10“® for J/0, (0.611 ^gi) ■ 10“® for T(IS'), (2.02 • 10“® for 

T(2S') and (2.44 ^}; 3 q) ■ 10“® for T(3S'). We find good agreement with the results reported by 
these authors except for the decay h T(IS') 7 , where their value is about a factor 7 smaller 
than ours. The reason is that we do not neglect the imaginary part of the direct contribution 
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to Ax(i 5 ) in fl42ll . which prevents |1 — Ax(i 5 )[ from becoming arbitrarily small. 

Our predictions may also be compared with the upper limits obtained from a recent first 
analysis of these rare decays reported by the ATLAS collaboration. They are Br(h —)■ J/0 7 ) < 
1.5 • 10-0 Br(h ^ T(l^) 7 ) < 1.3 • 10-0 Br(h ^ T(25) 7 ) < 1.9 • lO-^ and Br(h ^ T(3^) 7 ) < 
1.3 • 19 -^, all at 95% CL [^. It will require an improvement by a factor 500 to become 
sensitive to the h —)■ J/07 mode in the SM, while the SM branching fractions for the decays 
h —)■ T{nS) 7 are out of reach at the LHC. Nevertheless, as we will discuss below, these decay 
modes allow for very interesting new-physics searches. With 3 ab-^ of integrated luminosity, 
about 1.7 X 10® Higgs bosons per experiment will have been produced by the end of the high- 
luminosity LHC run Q. If the J/0 is reconstructed via its leptonic decays into muon pairs, 
the effective branching ratio in the SM is Br(h —)■ J/07 —)■ = 1.8 ■ 10-’^, meaning 

that about 30 events can be expected per experiment. If also the decays into e’''e- can be 
used, then ATLAS and CMS can hope to collect a combined sample of about 120 events. 
A detailed discussion of the experimental prospects for reconstructing these events over the 
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Figure 6 : Predictions (central values) for the ratios of the h V'y and h —)■ 77 
branching fractions with V = (f), J/ijj and T(IS') as functions of the rescaled Yukawa 
couplings normalized to the parameter k®® dehned in (HOI) . The black dots indicate the 
SM values. Coupling parameters inside the dashed white circle in the third plot are 
preferred by the current LHC data. See text for further details. 


background can be found in [10|]. Concerning the h —>■ 07 decay mode, a reconstruction 
efficiency = 0.75 was assumed for the 07 hnal state in [ll|, which appears to us as an 
optimistic assumption. In the SM one expects about dOOe^^y events per experiment in this 
mode, meaning that the two experiments can hope to look at a combined sample of several 
hundred events. Likewise, in the SM one expects about 2900epOy, events per experiment in the 
decay mode h —)■ 

In Figure[ 6 ]we show our predictions for the ratio of branching fractions (times 1000) dehned 
in flHTl) in the plane of the parameters ky/K®® and ky/K®®. We focus on the most interesting 
cases P = 0, J/0 and T(IS'). The corresponding plots for V = u would look very similar 
to that for P = 0 (apart from the overall scale of the branching fractions), while the plots for 
higher T(?t,S') resonances would look very similar to that for the T(IS') meson. For orientation, 
we mention that a value of 0.4 in these plots corresponds to a h —)■ P 7 branching fraction of 
about 10“®, assuming that the h —)■ 77 branching fraction is SM like. This assumption will be 
implicit whenever we quote absolute branching ratios below; otherwise the quoted numbers 
must be rescaled by Br(h —)■ 77 )/Br(h —)■ 77 )sm- The structure of our results flT2H implies that 
the rescaled Yukawa couplings always enter normalized to k®®. Hence, if a deviation from the 
SM is observed in any of these modes, then this could be caused either by a new-physics effect 
on the h —)■ 77 branching ratio (parameterized by k®®) or by a non-standard Yukawa coupling. 
The former effect would however be correlated among all decay channels. We observe that 
the h —)■ 07 branching ratio is rather insensitive to the CP-odd parameter k^, the reason 
being that this parameter enters quadratically and only via the direct contribution, which by 
itself is a small correction. An analogous statement holds (with less accuracy) for the case 
k ^ J/07. It is thus a reasonable approximation to study these decay modes under the 
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Figure 7: Predictions for the h ^ p'y and h ^ (p'y branching ratios, normalized to 
the h —)• 77 branching fraction, as functions of RpO ^ ( 2 /t„ + hrf)/3 and ~ Rg, 
respectively, normalized to The SM values are indicated by the red arrows. 


assumptions that = 0 and Kc = 0. The situation is different for the h —)■ T(nS ')7 decay 
modes, for which there is a strong cancellation between the direct and indirect contributions. 
The direct contributions are no longer a small correction, and hence the quadratic terms in ki, 
and Kb are important. The dashed white circle in the third plot indicates the current upper 
bound on the combination 



Kb 


2 

+ 


Kb 
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(46) 


To an excellent approximation, measures the deviation of the ratio Br(h —)■ bb) /Br(h 77 ) 
from its SM value. The Higgs bosons must be produced via the same production mechanism 
in both cases, so that possible new-physics effects in Higgs production cancel out. Since 
the h ^ bb mode is measured at the LHC in the rare Vh and tih associated-production 
channels, at present no accurate direct measurements of Xb-y are available. However, from 
the model-independent global analyses of Higgs couplings performed by ATLAS and CMS, in 
which all couplings to SM particles (including the effective couplings to photons and gluons) 
are rescaled by corresponding Ki parameters and also invisible Higgs decays are allowed, one 
obtains Xfyy = 0.63 ±0.27 for CMS (see Figure 17 in [^) and Xiyy = 0.67±0.32 for ATLAS (see 
Table 9 in j^). At 95% CL this (roughly) implies Xby < 1.3. Within this allowed region, the 
h —)■ T(IS') branching ratio varies by more than two orders of magnitude and can take values 
as large as 1.3 ■ 10“®. This might be accessible in the high-luminosity run of the LHC. If the 
T(IS') meson is reconstructed via its decays into muon or electron pairs, one could then hope 
for a sample of about 20 events with 3 ab“^ combining the ATLAS and CMS data sets. 

In order to better assess the theoretical uncertainties in our predictions, we now study 
the projections of the results onto the axis where the CP-odd couplings vanish. For the light 
mesons {V = p^,uj,(j)), setting ky = 0 has basically no impact on the branching ratios. In 
Figure [7] we show the ratio of branching fractions dehned in fl37)) as a function of the CP-even 
couplings Kv for h —)• p °7 and h —)■ ^ 7 . The width of the bands indicates the theoretical 
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Figure 8 : Predictions for the h and h —)■ T(IS') 7 branching ratios, normalized 

to the h —)■ 77 branching fraction, as functions of Kc and Kb, respectively, normalized 
to K^. The SM values are indicated by the red arrows. 


uncertainty. We have not included the small uncertainty in the values of the decay constants 
fv, because they can be eliminated using relation fl38|l . and we assume that by the time the 
h —)• 1/7 modes will be discovered the decay constants will have been measured more precisely 
than today. The corresponding plot for h wy would look identical to the left plot, but with 
a different vertical scale. While the theoretical uncertainties are small in all cases, we observe 
that the sensitivity of the branching ratios to the modihed Yukawa couplings is unfortunately 
rather week. For example, a hypothetical 10% measurement of the h —)■ branching ratio 

at the SM value would imply that \Kpo/Kf^\ < 0 . 8 , which is to say that a certain combination 
of the up- and down-quark Yukawa couplings is bounded not to exceed 80% of the 6 -quark 
Yukawa coupling. A 1% measurement would be required to obtain the more interesting bound 
\Kpo/K^\ < 0.08, which is still more than 100 times the SM value for RpO given in fl22D . The 
situation is not much better for the case h —>■ ^y. With a 10% measurement of the branching 
fraction at the SM rate, one would be able to conclude that —0.55 < R^/k^^ < 0.62. With 
a 1% measurement one would obtain the bounds —0.04 < R^/k^ < 0.08, which would come 
close to the SM value R^ ~ 0.02. Such a measurement is however out of the reach of the LHC. 

We now turn to the more interesting cases of radiative Higgs decays into heavy quarkonium 
states. In Figure [ 8 ] we show our predictions as a function of the physical parameters Kc (not 
Rc) and Kb, again assuming that the CP-odd couplings kc and kb vanish. In the latter case the 
impact of a possible CP-odd coupling on the branching fraction can be signihcant, and in the 
case of a measurement of a non-standard rate one should keep this possibility in mind. From 
the left plot in the hgure we conclude that a 20 % measurement of the h —)■ J/'ijj'y branching 
ratio at the SM value would allow one to constrain —0.51 < Kc/k^^ < 3.07, which would 
provide quite interesting information on the CP-even charm-quark Yukawa coupling. With 
a 10% measurement this range could be shrunk to 0.32 < k^/k"^ < 1.53, and with a 5% 
measurement one could reach 0.75 < Kc/k^^ < 1.19. Such accurate measurements serve as an 
interesting physics target for a future 100 TeV proton-proton collider. 
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Figure 9: Constraints on the effective coupling strengths Kfe/K®® and derived in 

two possible scenarios for future measurements of the ratios Br(h —)■ hh)/'Qi{h —)■ 77 ) 
(light blue) and Br(h —)■ T(IS') 7 )/Br(h —)■ 77 ) (orange). The allowed parameter space 
is given by the red-shaded intersection of the two rings. The black dot indicates the 
SM value. 


The situation with the h —)■ T {nS) 7 decay modes is different and quite interesting. In 
the SM the corresponding branching fractions shown in fH5|) are so small that these decays 
would be unobservable. The strong suppression arises from an almost perfect cancellation 
between the direct and indirect contributions to the decay amplitudes, which results from the 
fact that in the SM Re AT(n 5 ) ~ 1 within a few percent, see fH3|l . Thanks to this fortuitous 
fact, these decays offer a much enhanced sensitivity to the effects of new physics. For instance, 
the SM value of the h —)■ T(IS') 7 branching ratio of 4 x 10“® can be enhanced by a factor 
of more than 200 for Kb/K®® ~ —1 or Kb/zt®® ~ 3. The first of these possibilities would yield 
a h —)■ 66 rate consistent with current LHC measurements. For example, with a hypothetical 
25% measurement Br( 6 . —)■ T(IS') 7 )/Br( 6 , —)■ 77 ) = (0.4 ± 0.1) ■ 10“^ one would conclude 
from the figure that —1.21 < Kb/< —0.64, which would be a very significant piece of 
information and a spectacular sign of new physics. 

One may ask whether the current bounds obtained by the ATLAS collaboration already 
have a significant impact on the Higgs couplings. Unfortunately this is not the case. We find 
that the upper values reported in 2 l| imply approximately 
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(47) 


both at 95% CL. For comparison, we note that mt/rric ~ 268 and rrit/mb ~ 60. In other 
words, the current bounds derived from exclusive 6 . —)■ U 7 decays imply that the couplings of 
the charm and bottom quarks to the Higgs boson are not much stron ger than the top-quark 
Yukawa coupling (the more optimistic value K\ < 220 was quoted in (26|). 

We emphasize again that any experimental information on the rare radiative decays h —)■ 
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T{nS)'y should be interpreted in terms of an allowed region in the two-dimensional plane of 
the couplings and The one-dimensional projection shown in Figure | 8 ]is meant 

for illustrative purposes only. It is interesting to speculate about some possible scenarios that 
may arise at the end of the high-luminosity LHC run with an integrated luminosity of 3 ab“^. 
Existing estimates of the precision achievable on the Yukawa coupling to bottom quarks and 
the effective Higgs coupling to photons (our parameter |k®®|) suggest that, at 95% CL, the 
quantity Xiry dehned in flT 6 |) can be measured with a precision at least as good as 15% 1^. In 
Figure [9] we consider two possible future scenarios: 


(1) Ab.^ = 1.0 ±0.15, 

(11) Ab^ = 0.65 ±0.10, 


Br(h^ T(l^) 7 ) 
Br(h —)■ 77 ) 


< 0 . 2 - 10 -^ 


Br(h ^ T(lg) 7 ) 
Br(h —)■ 77 ) 


(0.4 ±0.2) ■ 10"% 


(48) 


In the hrst scenario the ratio Br(h —>■ bb)/BT{h 77 ) is measured at its SM value and an 
upper limit of about 0.5 • 10“® is placed on the h T(IS') 7 branching ratio. Even though no 
observation of this rare decay is accomplished, the assumed upper bound still has a non-trivial 
impact, as it limits the allowed values of Kb/n^ to the right half-plane. In particular, this 
would exclude the possibility that ~ — 1 - In the second scenario the parameter Ab^ is 

measured close to its current value but with higher accuracy, while a rough 50% measurement 
of the h —)■ T(IS') 7 branching ratio at about (0.91 ± 0.46) • 10“® is obtained. In this case one 
could exclude the SM value Kb/n^ = 1 and limit the allowed values of this ratio to the left 
half-plane. These speculative results nicely indicate the power of future searches for the rare 
exclusive decay h —)■ T(IS') 7 in the high-luminosity phase of the LHC. The statistics of such a 
search can be approximately doubled by including also the T(2S') 7 and T(3S') 7 hnal states. 

Several ideas for constraining the absolute value of the charm-quark Yukawa coupling in 
the high-luminosity run at the LHC or at a future 100 TeV collider have been put forward in 
[ioj ]. If such a measurement can indeed be made, it implies a circular allowed region centered 
at (0, 0) in the plane of the parameters and Kc/k®®. A measurement of the h — )■ J/'07 

branching ratio would be mainly sensitive to Kc/k"^ and hence conhne the couplings to a 
curved band, which intersects this region (see the center plot in Figure [ 6 ]) • In this way, it may 
be possible to perform an analysis similar to that shown in Figure [9] for the Higgs couplings 
to the charm quark. 

If at a future 100 TeV proton-proton collider one succeeds to collect large data samples 
of the rare decays h —)• V 7 —)■ with V = J/'ip oi T{nS), one might even speculate 

about the possibility to measure the polarization of both the vector meson and the photon by 
considering events in which the photon undergoes a nuclear conversion to an electron-positron 
pair, in an alog y with what was proposed for h —)■ 77 decay in and for the B —)■ K*'f 
process in (^. With such a measurement it would be possible to differentiate between the 
two structures in ([5]), which in 3-vector notation correspond to the products Sy-et^ and £y x e* 
of polarization vectors. In this way one would become sensitive to the sign of the ratio Kq/Kq 
of the CP-odd and CP-even Yukawa couplings, thus breaking the degeneracy beetween the 
upper and lower half-planes in Figure [9l 
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5 Conclusions 


We have performed a state-of-the-art analysis of the rare exclusive decays h —)■ I/7 in the 
context of generic extensions of the SM with modihed Higgs couplings. These decays are char¬ 
acterized by a destructive interference between two decay topologies. The direct contribution 
is governed by diagrams where the Higgs boson decays into a quark anti-quark pair, from 
which the vector meson is formed. It is proportional to the (modified) Yukawa cou plin g of 
the Higgs to that quark flavor. Using QCD factorization and techniques developed in |20[, we 
have derived a closed analytic expression for the direct amplitudes at next-to-leading order 
in as as an infinite sum over Gegenbauer moments renormalized at the scale /r ~ irth. In 
this way large logarithmic corrections are resummed, and the sensitivity of our predictions 
to hadronic input parameters is reduced. Power corrections to the factorization formula are 
suppressed by (AgcD/’^/i)^ or and can be safely neglected. The second, indirect 

contribution to the decay amplitude arises from diagrams where the Higgs boson decays into 
a photon and an off-shell neutral gauge boson, which converts into the vector meson. Due to 
the smallness of the Yukawa couplings, this indirect contribution is the dominant one in all 
cases but h —)• T(nS')7. In order to reduce the sensitivity to possible new-physics effects in 
the effective (loop-induced) ^77 vertex and thus get access to the quark Yukawa couplings, 
we consider the ratio of the h —?• U7 and h —?• 77 branching fractions, in which the indirect 
contribution drops out up to very small power corrections. In our analysis we account for the 
effects of flavor mixing, which can be important for the light mesons p°, uj and 0. 

We have derived numerical predictions for the h —)■ U7 branching fractions and studied 
the possibility of using such measurements as probes of new-physics effects on the light-quark 
Yukawa couplings. In the SM, the branching ratios we find typically lie in the range of few 
times 10“®, several orders of magnitude below the ATLAS bounds for the h —)■ J/'07 and 
h —)■ T(nS')7 branching ratios reported in j^. Nevertheless, the very high yield of Higgs 
bosons expected in the high-luminosity phase of LHC operation, combined with a dedicated 
experimental effort, could make measurements of these rare processes possible. We estimate 
that with 3ab“^ of integrated luminosity it will be possible to probe for 0{1) modifications 
of the real part of the charm-quark Yukawa coupling and 0(30) modifications of the real part 
of the strange-quark Yukawa coupling. 

We have emphasized that the decays h —)■ T (nS) 7 provide a golden opportunity to probe 
for new-physics effects on the bottom-quark Yukawa couplings. Due to a fortuitous cancellation 
between the direct and indirect contributions to the corresponding decay amplitudes, the SM 
branching fractions for these modes are more than two orders of magnitude smaller than the 
h —)■ J/'07 branching ratio. They are unobservable at the LHC and at any conceivable future 
collider. However, we show that with 0{1) modihcations of the 6-quark Yukawa couplings 
these branching fractions can be enhanced to an observable level. Any measurement of such 
a decay would be a clear signal of new physics. Moreover, as we have shown, a combined 
measurement of the two ratios Br(h T(nS') 7)/Br(h —?• 77) and Br(h —?■ 66)/Br(h —?■ 77) 
can provide complementary information on the real and imaginary parts of the 6-quark Yukawa 
coupling. This will allow one to probe the interesting option that k. —1 has the opposite 
sign as in the SM. 

Several extensions of our work are possible and worth pursuing. Decays into flavor off- 
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diagonal neutral mesons can be used to probe for flavor-violating couplings of the Higgs boson, 
which at tree level are forbidden in the SM. Another straightforward generalization is the 
application of our approach to the weak radiative decays h —?• M^Z and h —?• where 

M can be a pseudoscalar or vector meson, as well as decays into hnal states containing two 
mesons. This is left for future work. 

Exclusive radiative decays of the Higgs boson are an experimentally challenging endeavor, 
because the expected branching fractions are very small and the hnal states not easy to recon¬ 
struct. Nevertheless, we have argued that studies of these processes may not be impossible in 
the high-luminosity phase at the LHC and, even more so, at a future 100 TeV proton-proton 
collider. This would present us with a unique laboratory to study the real and imaginary 
parts of the Yukawa couplings of bottom and charm quarks, and probe for enhanced Yukawa 
couplings of the lighter quarks. We cannot think of any other way in which such direct studies 
of Yukawa couplings could be performed. 

Acknowledgments 

We are grateful to J. Brod, F. De Fazio, A. Grozin, U. Haisch, A. Kagan, G. Perez, F. Yu, 
J. Zupan and R. Zwicky for useful discussions. M.N. thanks the Institute for Theoretical 
Physics in Heidelberg for hospitality and support under a J. Hans D. Jensen Professorship. 
This work has been supported by the Advanced Grant EFT4LHG of the European Research 
Gouncil (ERG), the Gluster of Excellence Precision Physics, Fundamental Interactions and 
Structure of Matter (PRISMA - EXG 1098), grant 05H12UME of the German Federal Ministry 
for Education and Research (BMBF), and the DEG Graduate School Symmetry Breaking in 
Fundamental Interactions (GRK 1581). 


A Effects oi (jj — cj) mixing 


We express the physical mass eigenstates |a;) and |0) in terms of the flavor eigenstates |a;/) = 
[\uu) + \dd)) and | 0 /) = |ss) by means of the rotation 


Icn) 

10 ) 


/ cos 9^^ 
\ sin 9^^ 


- sin 9^^ 
cos 9^^ 



(A.l) 


In the limit where OZI-violating contributions are neglected, we can then relate the matrix 
elements of the flavor-specihc tensor currents in flT^ to decay constants dehned in terms of 
the analogous matrix elements 


10 ) = < - kvn , 

{Mk,£)\ |0) = (kV‘' - Ve*") 


(A.2) 


of the flavor eigenstates \ui) and 10/) with the corresponding flavor currents. Assuming isospin 
symmetry, this gives 


ru± _ rd± _ cos ^^0 _|_ 

J LJ J LJ \/2 ^ 


= -sin0^^/^^ 


(A.3) 
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and 


=COS 0 ^^/^^, 


rul. _ rdl. _ p± 

J d> J (b /TT J ujt ' 


'0 ''UJCp J (pj 1 J (j, J (p JuJi 

By definition, the qnantities and /f defined below flT^ are then given by 


= cos /j; + V2 sin , 

^_L ^ ^_L sin O^fp 
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From fll9p it then follows that 
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(A.7) 


In the limit of SU(3) flavor symmetry the ratio can be replaced by 1. Note that in 

general the mixing angle 9,^,p and the matrix elements in flA.2p may depend on the momentnm 
transfer If this is the case, the valnes of and entering in flA.Sp and flA.dp are different. 
All parameters entering the qnantities Sy in flA.Tp mnst then be evalnated at = my. 


B RG evolution equations 

The running quark masses and transverse decay constants obey the RG equation^ 
A ^ rugifi) = 7™(/x) rugifi) , p A = -7^(1/) /^(p), 


(B.l) 


where 7 ”^ and 7 ^ are the anomalous dimensions of the quark mass and the QCD tensor current 
in flT^ . At two-loop order these objects were hrst obtained in ^ and 6 l|, respectively. At 
NLO in RG-improved perturbation theory, the evolution equations have the solutions 


rugifi) = 

/v(a) = 


asifi) \ 2/3o 


as(Ao 
OisifJ'o] 


^ _ iTbo - a.(/x) - asifip) ^ 




dvr 


^g(Ao) 5 
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^0 ^ 


(B.2) 


^ ^ li bo - blip as(A) - «.(Ao) ^ 




dvr 
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0 ; 


^We follow the convention that the anomalous dimensions of coupling constants are defined with the opposite 
sign than those of operators. 
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The relevant one- and two-loop coefficients of the anomalons dimensions read 


7o"* = -QCf , 7™ = -3Cl - - CfCa + y , 

lo = , 7^ = —19Cp H—y CfCa —^ CfTfUj , 

where Cf = 4/3, Ca = 3 and Tp = 1/2 for SUc(3). Moreover, 


(B.3) 


l3o — -^C A 
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■TFUf, f^i = -c. 
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CaTf'^T'/ — ACfTf^^ 


(B.4) 


are the hrst two coefficients of the QCD /3-fnnction. Above, /xo is a low reference scale, at which 
the qnark masses and decay constants are calcnlated using some non-perturbative approach. 

The evolution of the leading-twist LCDAs at NLO order has been studied in joi-Ofil . 
Starting at two-loop order the scale dependence of aj/- (/i) receives contributions proportional 
to (/i) with k = 0,..., n. Dehning the vector a = (1, , • • •) containing the even 

Gegenbauer moments, one hnds a(/i) = f7(/r,/ro) where the evolution matrix U has a 

triangular structure with entries 


Unk (h) ho) 


■, k<n, 

^n^°(h,ho) ] k = n, 


(B.5) 


and Unki.k'i k-o) = 0 for A; > n (with even k,n > 0). Explicit expressions for E„(/r, po) and 
Po) can be found, e.g., in Appendix A of jo^. The relevant two-loop anomalous dimen¬ 
sions for the Gegenbauer moments of a transversely polarized vector meson were calculated in 


68 |, 1^. Numerically, we obtain for the evolution from the low reference scale po = 1 GeV up 


to the high scale /i = ruh 
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(B.6) 


where we have adjusted the number of light quark flavors when crossing the thresholds at 
H = rric and mf,. Due to the high value of the electroweak scale the mixing effects are strongly 
suppressed; for example, we obtain a^^{mh) ~ 0.304ap^(/io) — 0.007ay^(/ro) + 0.0008. When 
NLO evolution effects are included, the coefficients of the various terms shown in the second 
row of fl2^ get rescaled by factors of 1.001, 0.994, 0.984, 0.975, 0.969, respectively. Given 
that all present estimates of the hadronic parameters a)/- are afflicted with large theoretical 
uncertainties, it is sufficient for all practical purposes to use the leading-order solution fl25D . 












Meson V 

/#(2GeV)//v 

Method 

Reference 

P 

0.76 ±0.04 

lattice (unquenched) 

122] 


0.687 ±0.027 

lattice (unquenched) 

[71] 


0.742 ± 0.014 

lattice (quenched) 

[72] 


0.72 ±0.02 tofo 

lattice (quenched) 

[73] 


0.70 ±0.04 

light-cone sum rules 

[42] 


0.72 ±0.04 

our combination 


u 

0.707 ±0.046 

light-cone sum rules^ 

[42] 


0.71 ±0.05 

our combination 


0 

0.750 ± 0.008 

lattice (unquenched) 

Iffl 


0.780 ± 0.008 

lattice (quenched) 

[72] 


0.76 ±0.01 

lattice (quenched) 

[73] 


0.763 ± 0.041 

light-cone sum rules^ 

[42 


0.76 ± 0.04 

our combination 



Table 4: Compilation of theoretical predictions for the ratio fv at p = 2 GeV 

for light vector mesons. Values marked with a dagger are obtained by taking ratios 
of individual results for the two decay constants. In our combinations we adopt more 
conservative error estimates than in some of the original references. 


C Transverse vector-meson decay constants 


The direct contributions to the h ^ V'y form factors in fl2U]) involve the ratio /y(p-)//v of 
the vector-meson decay constants of tensor and vector currents, as dehned in ([HD and flT^ . It 
is reasonable to assume that non-perturbative evaluations of this ratio are subject to smaller 
theoretical uncertainties than calculations of the individual decay constants. For light mesons, 
predictions for the ratio of decay constants have been obtained using light-cone QCD sum rules 
and lattice QCD. Table 0] shows a compilation of available results, normalized at the reference 
scale /i = 2 GeV. The numbers in boldface show our own combinations of these results. 

For heavy quarkonia, the ratio of decay constants can be calculated using NRQCD. In¬ 


cluding the leading relativistic corrections [33| and one-loop QCD effects [32|, we obtain 


/v (h) 

fv 


mv 


2m, 
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1 - 3 {v^)v + 
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(C.l) 


where mg is the pole mass of the heavy quark, and the dots denote higher-order terms. 
Numerical values for the NRQCD matrix element {v‘^)v for fhe J/'ip and T(nS') states with 
n = 1,2,3 have been obtained from an analysis of the electromagnetic decays V —)■ e’''e“ 
including 0{as) corrections and the leading relativistic effects j^, The results obtained 
in this way, along with the adopted values of the pole masses of the heavy quarks, are compiled 
in Table O There are signihcant uncertainties related to the values of the heavy-quark pole 
masses and the NRQCD matrix elements which in our opinion have been underestimated 
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Meson R 

niQ [GeV] 

(uQv 

f^(2GeV)lfv 

J/t 6 

1.4 ±0.2 

0.225 IS:™ 

0.91 ±0.14 

T(1A) 

4.6 ±0.1 

-0.009 ±0.003 [7^ 

1.09 ±0.02 

T{2S) 

4.6 ±0.1 

0.090 ±0.011 [7^ 

1.08 ± 0.02 

T{3S) 

4.6 ±0.1 

0.155 ±0.018 [7^ 

1.07 ±0.03 


Table 5: NRQCD parameters for heavy quarkonia extracted from their electronic decay 
widths r(l/ ^ e+e“), and resulting values for the ratio of decay constants. 


in these analyses. It is well known that the concept of a pole mass is ill dehned beyond 


perturbation theorv and affected bv renormalon ambienities 76. 771. The values used in m 


are “one-loop pole masses”, but the intrinsic uncertainties in these values are parametrically 
of order Aqcd- Also, as emphasized in j^, the NRQCD expressions for the electromagnetic 
decay rates r(R —)■ receive very large perturbative corrections at two- and three-loop 

order, and this prevents an accurate extraction of the non-perturbative parameters {v^)v. In 
view of these remarks, the errors assigned on the 6 -quark mass and on the NRQCD matrix 
elements for the T{nS) states seem overly optimistic. In order to be conservative, we increase 
these errors by a factor of 2. This yields the results shown in Table [2j 


D Loop functions 


The loop functions describing fermionic and bosonic contributions to the off-shell h —)■ yR* 
decay amplitudes with R = 7 , Z have been derived hrst in ^|. In our notation, they reacfl 


Af{T,x) = 
Bf{T,x) = 


3r 


2x 


2(1 — x) 1 1 — X 

T 


9{^) -9{r/x) 


+ 1 


1 — X 


/(^) - 


1 — X L' 




2 x 


(d = 


1 — X 
1 

1 — X 


1 — X L' 


1 ~ -|- 


9{r)-9{r/x) 


+ 


3r 


3s^ I T 


(1 -x) 
2 x 


1 — X 


r 


(1 -x)^ 


2 (2 r) 2x(3 4s^) 


8x\ 

l 


/W - f{^/x) 


9{r) -9{x/x) 


(D.l) 


"'^We have corrected a typo in the expression for the function Af{T,x) given in eq. (4) of [51j . 
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where 


arcsin 


/(^) = < 


T 


I f, I + . 

-- In- , - tir 

4 V 


r>l, 

r < 1 , 


(D.2) 


and g{T) = r(l — r) /'(t). In the limit r —)■ oo one hnds that Af{T,x) —)■ 1, Bf{T,x) —)■ 1, 
x) ^7 and A'^ (r, x) -)> f - 7s^. 


E Coefficient functions Ay and Ay 


Here we list the complete expressions for the CP-even coefficients Ay and the CP-odd coeffi¬ 
cients Ay dehned in ([39]). For the former ones, we obtain 


Apo 


A,„ = 


A^ — 


Aj/p = 


[(0.068 ± 0.006) + i(0.011 ± 0.002)] + O.OOOl/ty. - O.OOOIk^ 




eff 

77 


[(0.068 ± 0.006) + i(0.011 ± 0.002)] - O.OOOIk^ - O.OOOlh^ 




eflf 

77 


[(0.093 ± 0.008) + z(0.015 ± 0.003)]+ 0m02Kw - 0.0002k, - O.OOOIk.^^ 


K 


eflf 

77 


[(0.281 ± 0.045) + i(0.040 ± 0.009)] it. 


K. 


eff 

77 


0.0004kvi/ — O.OOOOk., — O.OOOIk;, -|- O.OOOIk^ — O.OOOOk-^^ 


K, 


eff 

77 


A 


T(1S) 


[(0.948 ± 0.040) + i(0.130 ± 0.019)] k,, 


k! 


eff 

77 


0.019kvk — O.OOlKt — O.OOlzKr -|- (0.001 — 0.002 z)kc — O.OIOk-^z 

77ff ’ 

^77 

(E.l) 

where the parameters Ky for light mesons have been dehned in ffl^ . The corresponding 
expressions for the T(2S') and T(3S') mesons will not be shown explicitly, since they are very 
similar to that for the T(IS') state, see fl42p . In the above results the hrst term shows the 
direct contribution. The remaining terms, which originate from the power-suppressed h —)■ 
'-yZ* —)■ '-yV contribution and the effect of the off-shellness of the photon in the h —)■ 77 * jV 
contribution, are extremely small. Even for A^ and assuming a SM-like Higgs couplings to 
strange quarks, the theoretical uncertainty in the direct contribution is an order of magnitude 
larger than the power-suppressed terms. Only for T(IS') the power-suppressed terms reach 
the level of 10 “^, but still this contribution is smaller than the theoretical uncertainty in the 
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direct contribution. Our complete expressions for the CP-odd coefficients Ay are 

[(0.068 ± 0.006) i(0.011 ± 0.002)] + (0.429 0.003i)Kt (0.306 0.002i)K.^.^ 


Apo - rcp = 


..eff 

^77 


- rcp = 


(0.004 - 0.003i)K^ (0.002 - 0.002i)Kb -h (0.005 - 0.003 z)kc 

Zeff ’ 

Kyy 

[(0.068 ± 0.006) i(0.011 ± 0.002)] + (0.429 0.003z)Kt (0.306 -h 0.002 z)k. 


'77 


^efF 


A<a - rcp = 


(0.004 - 0.003i)K^ (0.002 - 0.002i)Kfe -h (0.005 - 0.003i)iCc 

Tiff ’ 

^77 

[(0.093 ± 0.008) -h i(0.015 ± 0.003)] k^ + (0.429 -h 0.003i)Kt -h (0.306 -h 0.002i)At. 


'77 


.eflf 

''77 


Aj/^ - rcp = 


(0.004 - 0.003i)K^ (0.002 - 0.002i)Kb (0.005 - 0.003 z)kc 

Tiff ’ 

^77 

[(0.277 ± 0.045) i(0.043 ± 0.009)] Kc (0.429 0.003z)Kt -h (0.306 0.002i)K, 


'77 


A 


T(15) - ''’CP 


^eff 

(0.004 - 0.003i)K^ (0.003 - 0.002i)Kb 

Teff ’ 

^77 

[(0.945 ± 0.040) i(0.132 ± 0.019)] kb -f (0.427 -h 0.003i)icf -h (0.306 -h 0.002i)it, 


'77 


^eff 

77 


(0.004 — 0.002i)KT- -|- (0.004 — 0.001i)Kc + 0.010/l.yZ 


^eff 

(E.2) 

It is a good approximation to only keep the direct contributions, which are likely to give rise 
to the dominant effects. 
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